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ABSTRACT 

The uninterrupted 7-day ASCA observations of the TeV blazar Mrk 421 in 1998 have clearly revealed 
that X-ray flares occur repeatedly. In this paper, we present the results of the time-resolved spectral 
analysis of the combined data taken by ASCA, RXTE, BeppoSAX, and EUVE. In this object - and in 
many other TeV blazars - the precise measurement of the shape of the X-ray spectrum, which reflects 
the high energy portion of the synchrotron component, is crucial in determining the high energy cutoff of 
the accelerated electrons in the jet. Thanks to the simultaneous broadband coverage, we measured the 
0.1-25 keV spectrum resolved on time scales as short as several hours, providing a great opportunity to 
investigate the detailed spectral evolution at the flares. By analyzing the time subdivided observations, 
we parameterize the evolution of the synchrotron peak, where the radiation power dominates, by fitting 
the combined spectra with a quadratic form (where the vF v flux at the energy E obeys \ogvF v {E) = 
log(i/F UtPca k) — const x (logE — log Speak) 2 ). In this case, we show that there is an overall trend that the 
peak energy -E pca k and peak flux vF v ^ ea \^ both increase or decrease together. The relation of the two 
parameters is best described as -E p0 ak oc vF®-£ eak for the 1998 campaign. Similar results were derived 
for the 1997 observation, while the relation gave a smaller index when included both 1997 and 1998 
data. On the other hand, we show that this relation, and also the detailed spectral variations, differ 
from flare to flare within the 1998 campaign. We suggest that the observed features are consistent with 
the idea that flares are due to a appearance of a new spectral component. With the availability of the 
simultaneous TeV data, we also show that there exists a clear correlation between the synchrotron peak 
flux and the TeV flux. 

Subject headings: BL Lacertae objects: individual (Mrk 421) — galaxies: active — radiation 
mechanisms: non-thermal — X-rays: galaxies 



1. INTRODUCTION 

The uninterrupted 7-day ASCA observation of the TeV 
blazar Mrk 421 in 1998 has revealed that day-scale X-ray 
flares seen in previous observations were probably unre- 
solved superpositions of many smaller flares (Takahashi et 
al. 2000). The nearly-continuous observation allowed not 
only the possibility to track the individual flares entirely 
from the rise to decay, but it also enabled quantitative 
statistical tests of the time series by employing the power 
spectrum or the structure function (Kataoka et al. 2001; 
Tanihata et al. 2001). 

The main characteristic of blazars is their high flux ob- 
served from radio to 7-rays coupled with strong variability 
and strong polarization. These properties are now suc- 
cessfully explained by the scenario where blazars are ac- 
tive galactic nuclei (AGN), possessing jets aligned close 
to the line of sight, and accordingly the Doppler-boosted 
non-thermal emission from the jet dominates other emis- 
sion components (e.g. Blandford & Konigl 1979; Urry & 
Padovani 1995). This is what makes blazars critical in 
understanding jets in AGN. 

The broadband spectra of blazars consist of two peaks, 
one in the radio to optical-UV range (and in some cases, 
reaching to the X-ray band), and the other in the hard 
X-ray to 7-ray region. The high polarization of the radio 



to optical emission suggests that the lower energy peak is 
produced via the synchrotron process by relativistic elec- 
trons in the jet. The higher energy peak is believed to 
be due to Compton up-scattering of seed photons by the 
same population of relativistic electrons. Several possibil- 
ities exist for the source of the seed photons; these can be 
the synchrotron photons internal to the jet (Jones, O'Dell, 
& Stein 1974; Ghisellini & Maraschi 1989), but also exter- 
nal, such as from the broad emission line clouds (Sikora, 
Begelman, & Rees 1994) or from the accretion disk (Der- 
mer, Schlickeiser, & Mastichiadis 1992; Dermer & Schlick- 
eiser 1993). 

The blazars with peak synchrotron output in the X-ray 
range also emit strongly in the 7-ray energies, and the 
brightest of those have been detected in the TeV range 
with ground-based Cherenkov arrays. These are the so- 
called "TeV blazars." In TeV blazars, variability of the 
synchrotron flux is measured to be the strongest and most 
rapid in the X-ray band, and thus it provides the best 
opportunity to study the electrons that are accelerated to 
the highest energies. In particular, the synchrotron peak is 
a very important observable in two aspects: first because 
the flux at the peak represents the total emitted power 
from the blazar, and secondly because the peak frequency 
reflects the maximum energy of radiating particles gained 
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in the acceleration process. 

Mrk 421 is among the closest known blazars, at red- 
shift of 0.031, and was the first blazar (and also the first 
extragalactic source) discovered to be a TeV emitter. It 
was first detected as a weak source by EGRET (Lin et al. 
1992), and 9 months later, Whipple detected a clear sig- 
nal from this object between 0.5 and 1.5 TeV (Punch et 
al. 1992; Petry et al. 1996). Flux variability on various 
time scales has been observed, including a very short flare 
with a duration of ~1 hour (Gaidos et al. 1996). Ever 
since, Mrk 421 has been repeatedly confirmed to be a TeV 
source by various ground-based telescopes. It has also been 
one of the most studied blazars, and have been a target of 
several multi-wavelength campaigns (Macomb et al. 1995; 
Maraschi et al. 1999; Takahashi et al. 2000). 

The multi- wavelength campaign of Mrk 421 in 1998 
was one of the first opportunities to observe a blazar 
in the TeV range using several telescopes located in dif- 
ferent locations in the world, so as to have as continu- 
ous coverage as possible (Aharonian et al. 1999; Piron 
et al. 2001). Observations in other frequencies included 
X-ray observations by ASCA, RXTE (Takahashi et al. 
2000), and BeppoSAX (Maraschi et al. 1999; Fossati et al. 
2000a, b), EUV observations by EUVE, optical observation 
with BVRI filters organized by the WEBT collaboration 
(http://www.to.astro.it/blazars/webt/), and radio obser- 
vations at the Metsahovi Radio Observatory. 

In this paper, we present the results of the spectral anal- 
ysis of X-ray and EUV data during the 1998 April cam- 
paign. In particular, we estimate the location of the syn- 
chrotron peak in the vF v spectrum. So far, quantitative 
studies of the variation of the synchrotron peak have been 
conducted only for the two brightest blazars, Mrk 421 and 
Mrk 501. The largest variation of the synchrotron peak 
energy was observed in Mrk 501 by BeppoSAX, where the 
peak energy shifted ~2 orders of magnitude (Pian et al. 
1998). Collecting data from different epochs (separated 
by as long as years), Tavecchio ct al. (2001) have shown 
the relation of the form -Ep ea k °c -Fo\_iookcV w ^ n 71 — ^ 
A similar analysis was done by Fossati et al. (2000b) for 
the Mrk 421 data obtained by BeppoSAX in 1997 and 
1998, which showed a relation of i? pC ak oc ^o\_iok e v wrtn 
n=0.55. In this paper we describe a continuous 7-day vari- 
ation of the synchrotron peak, which allows us to inves- 
tigate the dynamical change of the synchrotron spectrum 
during the flares. 

We first describe the observations and data analysis in 
§2. The results of spectral analysis are described in §3; 
a summary of the results and a discussion of our findings 
with the emphasis on the temporal evolution of the spec- 
trum are presented in §4. 

2. OBSERVATION AND DATA ANALYSIS 

The week- long multi- frequency campaign of Mrk 421 
was carried out in April 1998 (Takahashi et al. 2000). The 
goal of this paper is to study the evolution of the high 
energy synchrotron spectrum as a function of time dur- 
ing the campaign, and also the correlation with the TeV 
flux. This requires the best possible spectral coverage in 
the X-ray range, for which we assembled all available data 
collected during the campaign, including: pointings by X- 
ray observatories ASCA, RXTE, and BeppoSAX, as well 



as the EUV observations by EUVE. 

2.1. ASCA 

Mrk 421 was observed with ASCA for 7 days continu- 
ously (except for the Earth occultation) during 1998 April 
23.97 - 30.8 UT (PI: T. Takahashi). The ASCA detectors 
- two SISs (Solid-state Imaging Spectrometers; Burke et 
al. 1991; Yamashita et al. 1997) and two GISs (Gas Imag- 
ing Spectrometers; Ohashi et al. 1996) were in operation. 
We refer to Tanihata ct al. (2001) to the details of the ob- 
servation modes and data reduction. The obtained spectra 
were rebinned such that all bins have the same statistics. 

Since late 1994, the efficiency of both SIS detectors be- 
low ~1 keV has been decreasing over time. The problem 
is believed to be caused by the increased residual dark cur- 
rent level, and also by the decrease in the charge transfer 
efficiency, although effort is still underway to model the 
effects. While it has been reported that it may be possible 
to parameterize the degradation by an additional absorp- 
tion column as a function of time (Yaqoob 2000), since our 
observation concerns the precise shape of the continuum 
spectrum, we will use only the data obtained by the GIS 
for spectral analysis. (Note that the X-ray spectra in Fig- 
ure 2 of Takahashi et al. (2000) show the data obtained by 
SIS. For this, there is a too strong decrease of inferred flux 
at the lower energy, which is due to instrumental effects 
described above.) 

2.2. RXTE 

The RXTE observations (PI: G. Madejski) were mainly 
coordinated to coincide with the ASCA observation pe- 
riod, but were conducted mostly only during the low- 
background orbits. The RXTE extends the ASCA band- 
pass to a higher energy with the instruments PCA (Pro- 
portional Counter Array; 2-60 kcV) and HEXTE (High 
Energy X-ray Timing Experiment; 15-200 keV). However, 
because both the flux from the source and the sensitivities 
of the detectors decrease towards higher energies, in the 
analysis below we use only the data from the PCA, in the 
energy range of 3-25 keV. 

All data reduction was performed using the HEASOFT 
software packages, based on the REX script, provided by 
the RXTE Guest Observer Facility (GOF). The screen- 
ing criteria excluded the data when the elevation angle 
from the Earth's limb was lower than 10°, and the data 
acquired during the SAA passages. Since the PCUs are 
activated during the SAA passage, we used the data only 
after when the background level dropped to the quiescent 
level and became stable, which is typically 30 minutes after 
the peak of a SAA passage. 

For the best signal-to-noise, we selected events only from 
the top xenon layer, and excluded events from other gas 
layers. For the maximum consistency among the observa- 
tions, we used only the three detectors which were turned 
on for all the observations (known as PCU0, PCU1, and 
PCU2). Since the PCA is not an imaging detector, the 
blank-sky observations must be used to estimate the back- 
ground spectrum. For this, we used the background mod- 
els provided by the RXTE GOF to evaluate the back- 
ground during each observation. We note that the residual 
background contamination of the data within the energy 
bandpass selected by us was always significantly lower than 
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the source counts. The mean net RXTE PCA count rate 
from Mrk 421 in the 10 - 20 keV range for 3 PCUs as 
above was at least <~ 3 counts s" 1 while the total residuals 
from blank-sky in this energy range plus the fluctuations of 
the Cosmic X-ray Background were lower than 0.2 counts 
s- 1 (see HEASARC Web page regarding the RXTE PCA 
background), providing the confidence that an imprecise 
background subtraction did not affect our results. 

2.3. BeppoSAX 

The BeppoSAX observations (PI: L. Chiappetti) started 
3 days before the ASCA observations (Maraschi et al. 1999; 
Fossati et al. 2000a, b). We use the data obtained by the 
LECS (Low Energy Concentrator Spectrometer; 0.1-10 
keV), which extends the bandpass to lower energy, and 
also the MECS (Medium Energy Concentrator Spectrom- 
eter; 1.3-10 keV). The data analysis was based on the 
linearized event files, together with the appropriate back- 
ground event files and response matrices, provided by the 
BeppoSAX Science Data Center (SDC). Data reduction 
was performed using the standard HEASOFT software, 
XSELECT. 

We extracted the source photons from a circular region 
centered at the source with a radius of 8 arcmin for both 
LECS and MECS. For the background photons, we con- 
structed the spectra by extracting the photons from the 
same detector regions using blank-sky observations. Each 
spectrum was rebinned using the grouping templates pro- 
vided by the BeppoSAX SDC, and we limited the energy 
range to 0.12-4.0 keV for LECS, and 1.65-10.5 keV for 
MECS, as suggested in the Cookbook for BeppoSAX NFI 
Spectral Analysis (1999). 

2.4. EUVE 

The EUVE observations covered the entire campaign of 
Mrk 421 from April 19.4 to May 1.2 (PI: T. Takahashi). 
The source photons were extracted in a 12 arcmin aper- 
ture. In order to estimate the absolute flux in the EUV 
range, we normalized it to be consistent with the simul- 
taneous data obtained by the other observatories. Given 
that the EUV spectrum is difficult to determine within the 
rather limited bandpass, we used only the observed inte- 
grated count rate (cf. Marshall, Fruscione & Carone 1995). 
For this, we used the period where the observations by all 
observatories overlapped (MJD 50927.03 - 50927.34; §3.1). 
During this period, the average EUVE count rate was 
0.46±0.01 cts/s. At the same time, the 0.13-0.18 keV flux 
derived from the combined LECS, MECS, CIS, and PCA 
spectrum gives vF v = (1.78 ±0.04) x 10~ 10 erg cm~ 2 s _1 , 
where the error represents the statistical error. 

An additional uncertainty arises from the fact that the 
difference in the spectral shape is not taken in to account. 
This can be estimated from the results in Marshall, Fr- 
uscione & Carone (1995). Taking A H = 1.5 x 10 20 cm" 2 
for Mrk 421, assuming that the spectrum index varies from 
0.5 to 1.0, the flux density would vary from 1363-1394 /zJy 
for 1 cts/s in the EUVE range. This gives an additional 
2.5% error as a systematic error. 

In the following analysis, we use the above value for 
scaling the observed EUVE count rate to flux. Thus, 

vF„(EUVE) = 1.78 x 10~ 10 ergs^^m" 2 , (1) 



where F cts / S is the EUVE count rate. 

3. RESULTS 

3.1. The Shape of the 0.1-25 keV Spectrum 

During April 24, observations by all 3 X-ray satellites, 
ASCA, BeppoSAX, RXTE, and aiso EUVE happened to 
exactly overlap continuously for 26.5 ksec (MJD 50927.03 
- 50927.34). This provides us with high quality data 
over a very broad energy coverage, spanning more than 
2 decades, from 0.12 to 25 keV. The net exposure for each 
instrument during the overlapping 26.5 ks period were 9.1, 
13.3, 3.6, and 5.5 ksec for CIS, PCA, LECS, MECS, re- 
spectively. 

We first fitted the spectrum for each instrument with a 
power-law model, with a low energy absorption fixed to 
the Galactic value, iV H = 1.5 x 10 20 cm" 2 (Elvis, Wilkes, 
& Lockman 1989). The results are shown in the top half of 
Table 1, suggesting that the fit is not adequate for any of 
the individual instruments. This certainly suggests that a 
more complex model, rather than just a simple power-law 
function, is required to describe the spectrum. 

The shape of the synchrotron spectrum is determined 
from the energy distribution of the radiating electrons, 
which cuts off at the high energy representing the accel- 
eration limit. As described in the introduction, the high 
energy end of the synchrotron spectrum in TeV blazars is 
located in the X-ray range; this is most likely the reason 
that the high-quality X-ray data cannot be described by 
a single power-law. 

Given the curved spectrum for each instrument, we fit- 
ted the combined X-ray spectrum with various functions 
to model the shape of the observed curvature. We note 
that our aim is not to develop the shape of the curva- 
ture from first principles, but rather, to reproduce the 
observed data. For this reason, we considered the follow- 
ing models for comparison: (1) A broken power-law, (2) 
a curved power-law, (3) a power-law with an exponential 
cutoff, and (4) a quadratic function in log v-\ogvF v plane. 
The curved power-law that we define below is similar to 
the broken power-law model, but the change in the spec- 
tral index is not discrete. The form attempted by us is 
F(E) = KE- Fl (l + t^)- V2 - The form for the exponen- 
tially cutoff power law is F(E) = KE- Tl exp(-E/E hl ), 
and the quadratic function is described as logv F V (E) = 
log(yF VtPeek ) - a(log£ - log.Bpoak) 2 . 

Since there are slight differences in the cross-calibration 
between the instruments, we allow the normalization to 
vary for each instrument. For all results below, the flux 
will be normalized to the GIS value, unless otherwise 
noted. The iVn is fixed to the Galactic value. 

By fitting the spectrum with the models above, we first 
found that the power-law with an exponential cutoff is not 
a good representation of the data: the model cutoff is too 
sharp at high energies. The broken power-law model is also 
not acceptable, and the curved power-law is significantly 
better. We also find that the quadratic function fits the 
spectrum surprisingly well. The fitted spectrum and the 
residuals to the quadratic model are shown in Figure 1, 
where all LECS, MECS, GIS, and PCA spectra are shown 
to be well described with this single function. The results 
of the spectral fits by the 4 models are summarized in in 
the bottom half of Table 1. 
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3.2. The Time Resolved Spectra 

The EUVE observation covered the whole campaign 
(April 19.4 to May 1.2), and thus completely overlapped 
our 7-day ASCA observation. This adds a data point 
at 0.13-0.18 keV in addition to the ASCA spectrum in 
the 0.7-7 keV range, providing us the best opportunity 
to study the continuous spectral variability. In fact this 
is the longest continuous monitoring of this object so far 
spanning over the 0.13-7 keV range. 

We divided the whole observation into shorter segments 
with each having a duration of 10 ksec. In order to com- 
bine the GIS and EUVE data, the GIS count rate spectra 
must first be converted to the source flux spectra as the 
EUVE count rate data converts directly into flux (§2.4). 
For this, we first fit the GIS spectrum with a curved power- 
law model with Galactic absorption. This gives the source 
spectrum in the GIS energy range. Note that this pro- 
cedure depends on the selected model, but long as it de- 
scribes the spectrum well, this should introduce the small- 
est bias. We then combine this with the EUVE data, 
which gives the combined source spectrum. 

In fitting the combined EUVE /GIS spectra, we first 
started with a single power-law model. In similarity to the 
results in §3.1, none of the 56 spectra could be fit with a 
single power-law function, with their average \t of 7.7, im- 
plying the curvature of the spectra for all data sets. Given 
that a quadratic model fits the spectrum equally well as 
the curved power-law model, but with fewer parameters 
(see §3.1), we will assume the quadratic model for all data 
sets hereafter. In fact, because of the energy gap between 
the EUVE and GIS, the fit with a curved power-law model 
resulted in having multiple minima in the x 2 plane (since 
the parameters are correlated) , and the exact value of the 
best-fit parameters depended highly on the initial value. 

We checked whether the fit was adequate for all data seg- 
ments. The xl f° r eacri n t (shown in the bottom panel in 
Figure 2) ranged from 0.71 to 1.45. The number of degrees 
of freedom ranges from 60 to 187, which gave P(x 2 ) >5% 
for 53 out of the total of 56 data points. For the other 
3 data points (corresponding to times 230-240, 270-280, 
and 450-460 ksec) , we found several line or edge- like struc- 
tures in the spectra, which resulted in the increase of xt- 
Although it is not clear whether this is due to some system- 
atic effect or if it is a true property of the source, since our 
aim is to model the continuum spectrum correctly, we ex- 
cluded the energy ranges where those features are present. 
This improved the fit significantly, which is shown as the 
crosses. With this, all the 56 spectra were well described 
with the quadratic function. 

The time evolution of the derived parameters in the fit 
is shown in the top 3 panels in Figure 2; the peak flux, 
peak energy, and the curvature parameter a (see defini- 
tion in Table 1). The peak energy is observed to shift 
between 0.5 keV and 1.2 keV, in which a general trend is 
found such that the peak energy is higher when the peak 
flux is higher. This is demonstrated as the filled circles in 
Figure 3. The fit results are listed in Table 2. 

We do note that the actual validity of using a quadratic 
function in the entire 0.1-25 keV range was checked in de- 
tail (including an examination of trends in residuals) only 
for the epoch when the overlapping BeppoSAX, ASCA, 
and RXTE observations were available. Since we mea- 



sured only the EUVE flux below 0.7 keV in the combined 
EUVE/GIS fit, the derived peak location can depend on 
the model function we assume. For instance if any function 
was allowed, the peak can be located at any energy within 
the gap. For this, we emphasize that the derived peak lo- 
cation is the estimated peak, assuming that the spectrum 
gap can be extrapolated with the same quadratic function, 
and also that the error bar describes only the statistical 
error. 

We also analyzed the BeppoSAX data during the cam- 
paign. Furthermore, in order to make a comparison with 
data obtained from another campaign, we also analyzed 
the BeppoSAX data obtained in 1997 (PI: G. Vacanti, 
L. Chiappetti). We divided the observation into seg- 
ments of the same 10 ks duration, and fit the combined 
LECS/MECS spectra with a quadratic function. For the 
BeppoSAX data, we also fit the spectra with a curved 
power-law model for comparison. The derived peak flux 
and peak energy from both models are listed in Table 3. 
In fitting the combined LECS/MECS spectra, the normal- 
izations are allowed to vary independently, and the values 
in Table 3 are normalized to the MECS2 value. The fitting 
results show that although the spectral fit parameter is one 
less in the quadratic function, both models equally repre- 
sent the observed spectra. In contrast to the combined 
EUVE/GIS data, there is no energy gap in the BeppoSAX 
data, which provides support that our assumed quadratic 
form is appropriate in estimating the peak location in the 
EUVE/GIS case. The derived peak flux and peak energy 
using both functions are similar, but we also note that the 
quadratic function gives a slightly lower peak energy than 
the curved power-law function. 

The derived peak flux and peak energy assuming the 
quadratic model for the BeppoSAX observations are shown 
in Figure 3 as the open circles (for 1998) and squares (for 
1997). Note that the BeppoSAX data points in Figure 3 
are multiplied by a normalization factor of 1.07, which we 
adopted from the combined fit in §3.1 It can be seen that 
the results during the 1998 campaign lie on a similar slope 
as the results derived from the EUVE/GIS data during 
the same campaign. 

The relation of peak flux vF v>pea k and peak energy 
Epeak for the EUVE/GIS data set can be best described 
with a power-law function of E poa k oc f / E" pcak , where 
a=0.72±0.02 (or e=0.76±0.02 in form E pcak oc E £ A _ wkeV ; 
see description below). By including also the BeppoSAX 
data points during 1998, the index becomes a=0.77±0.02 
(e=0.72±0.02). It can be seen that the data from 1997 lie 
on a somewhat different slope. The relation derived from 
the 1997 data alone gives a=0.96±0.09 (e=0.79±0.06), 
slightly steeper than the index derived from the 1998 cam- 
paign. Fitting all data together we obtain a=0.59±0.01 
(e=0.50±0.01). All errors correspond to a la error. 

A similar analysis was performed by Fossati et al. 
(2000b) for the 1997 and 1998 BeppoSAX data. Fossati et 
al. (2000b) note a tight relationship between the measured 
integrated 0.1-10 keV flux Eo.i_iokcV arid the peak energy 
Epeak such that Ep C ak oc Fo.i-iokeV witn £=0.55±0.05, 
similar to the index derived from our results including 
all data sets. For a direct comparison, we fit simultane- 
ously both 1997 and 1998 for BeppoSAX data points only. 
This gave a=0.52±0.02 (e=0.45±0.01). The estimated e 
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is slightly smaller than reported in Fossati et al. (2000b), 
which might be affected by different subdividing of the ob- 
servations. Nonetheless, it is worth noting that the slope 
derived from the combined data appears be softer (smaller 
a or e) than each epoch considered alone. 

Furthermore, the relation appeared to also vary between 
individual flares during the 1998 campaign. In some flares 
the peak energy clearly increases together with flux, while 
the correlation is weaker in others. This can be seen clearly 
in Figure 4, which shows the evolution of the synchrotron 
peak in each flare. The peak energy has a trend to shift 
to the higher energy and then return to the initial value 
together with the flux, but there is no common slope. 

3.3. Spectral Variations at Flares 

In the previous section, we described the spectral vari- 
ations as the evolution of the location of the peak in the 
vF v spectrum. Here, we take advantage of the wide energy 
range coverage afforded by the use of by EUVE, ASCA, 
and RXTE data to investigate spectral variations in more 
detail. Given that the RXTE observations happen mostly 
during the flare minima rather than at their peaks (see 
Figure 1 in Takahashi et al. 2000), here we focus on the 
beginning (rising portion) of the flares. 

We divided the observations into 10 ksec segments, sim- 
ilar to the previous section. For each segment, the GIS 
and PCA spectra are fitted simultaneously with a curved 
power-law function with an additional constant factor for 
the normalization. The derived source spectra in the GIS 
and PCA range are combined with the EUVE data point, 
which gives the total source spectrum in the energy range 
of ~0.1-25 kcV. In the top panels of Figures 5-7, we show 
the evolution of the uF v spectrum derived in this manner 
for the 2nd, 4th, and 5th flare (starting from 80 ksec, 260 
ksec, and 340 ksec, respectively, as defined in Figure 2). 
Note that the spectra are further rebinned for this Figure. 
For each flare, we also plot a subsequent spectrum, but 
those only include the EUVE and GIS data. 

To see the evolution more clearly, we calculated the ra- 
tio of the rebinned spectra by dividing these "subsequent" 
spectra by the spectrum measured at the onset of a given 
flare (integrated over 80-90, 260-270, 340-350 ksec, re- 
spectively). This is shown in the bottom panels of Fig- 
ures 5-7. Interestingly, these three flares all showed dif- 
ferent features. It can be seen that Flare 2 starts with a 
hardening at the very highest energy range with a concave 
shape in the ratio. As time progresses, the whole energy 
range becomes dominated by a new spectral slope. In con- 
trast to this, in Flare 4 the lower energy photons appear 
to precede the total flux rise. Flare 5 was similar to Flare 
2, but all energy bands increase together, with a larger 
amplitude at higher energies. We also studied the spec- 
tral variation for other flares, but a detailed investigation 
was rather difficult without a denser temporal coverage by 
RXTE. 

3.4. The X-ray and TeV Correlation 

Since X-ray and TeV energies both reflect the electrons 
accelerated to the highest energies in TeV blazars, simul- 
taneous X-ray and TeV observations provide an important 
tool to probe the emission region of the jet (e.g., Krawczyn- 
ski et al. 2002, and references therein). On the other hand, 



since blazars are highly variable on short time scales, ex- 
actly simultaneous observations are essential. Thanks to 
the relatively high brightness of the source during the 1998 
campaign, it was possible to generate a significant mea- 
surement for the separate TeV pointings every night (Aha- 
ronian et al. 1999; Takahashi et al. 2000). Here we provide 
the X-ray spectrum from the exactly simultaneous time 
regions with the TeV observations. 

For each TeV observation by HEGRA (Aharonian et al. 
1999) and Whipple (Takahashi et al. 2000), we collected 
all existing X-ray data during the TeV exposure, from the 
ASCA, RXTE, and EUVE observations (the CAT obser- 
vations mostly overlap with the HEGRA observations). In 
particular, all Whipple pointings during the ASCA obser- 
vations happened to coincide with the RXTE observations, 
providing an X-ray spectrum ranging from 0.1-25 keV, si- 
multaneous with the TeV observation. 

We generated the combined X-ray spectrum using the 
same procedures as described in the previous sections. We 
then fitted the X-ray spectrum for each period assum- 
ing a quadratic function. The X-ray spectra during the 
HEGRA observations (consisting of the EUVE and ASCA 
data) are shown in Figure 8(a). The spectrum for each ob- 
servation is shown as different symbols, and the dashed line 
shows the modeled quadratic function derived via the fit. 
The X-ray spectra during the Whipple observations, which 
consist of the EUVE, ASCA and RXTE data are shown 
in Figure 8(b). Note that the spectra are further rebinned 
in the figures. While each integration time is as short as 
several hours, it can be seen that the curved shape of the 
high energy peak of the synchrotron spectrum is clearly 
resolved, and it is evident that the peak is located in the 
X-ray range. 

The results of fits to each X-ray spectrum concurrent 
with the TeV pointings are summarized in Table 4. We 
also list the simultaneous TeV fluxes, in which the HEGRA 
data were taken by Aharonian et al. (1999). In Figure 9, 
we show the correlation of the derived location of the syn- 
chrotron peak flux with the TeV flux. Since the sensitive 
energy ranges of the two TeV telescopes are slightly dif- 
ferent, we used the integrated flux in Crab units. A good 
correlation can be seen, showing that the TeV flux is higher 
when the synchrotron peak flux is higher. A best fit to the 
relation with a power-law function gives Ft v oc F s " poak 
with a=1.7±0.3. 

4. SUMMARY AND DISCUSSION 

We have performed detailed analysis of the combined 
EUVE, ASCA, BeppoSAX, and RXTE data collected dur- 
ing the long look campaign of Mrk 421 in 1998, resulting in 
the measurement of the time-resolved spectrum in a broad 
energy range of 0.1-25 keV in segments as short as several 
ksec. These are among the highest quality spectra with 
regard to both photon statistics and energy coverage so 
far for any blazar, providing the precise spectral shape at 
the highest energy end of the synchrotron spectrum. 

We have shown that both curved power-law function 
and quadratic function in log v-\o i gvF v space reproduce 
the combined EUVE, BeppoSAX, ASCA, and RXTE spec- 
trum. For the epochs where BeppoSAX data do not exist, 
we assumed that the spectral gap between EUVE and 
ASCA can be extrapolated with a quadratic function, and 
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we showed that the energy and the flux at the synchrotron 
peak in the vF v spectrum are correlated, showing an over- 
all trend of a higher synchrotron peak energy for higher 
peak flux, but the details of this relation differed from flare 
to flare. In particular, within the three flare beginnings, 
one flare started from a rise in the higher energy while an- 
other started from a rise in the lower energy. The relative 
amplitude of the rise in different energies also differed from 
flare to flare. An interesting result was shown from the 2nd 
flare, where it was indicated that this flare started with a 
hardening at the higher energy. This feature is difficult 
to describe by a simple cooling or acceleration of a sin- 
gle electron distribution (see e.g. Kataoka 2000), and thus 
this concave curvature indicates that two different spectral 
components were likely to co-exist at the beginning of the 
flare. We remark that such a behavior is observed in only 
one flare, but the importance is that we observed a flare 
that could not be explained by a single electron distribu- 
tion. This suggests an appearance of a new component, 
which generated the observed flare. 

Recently, many theoretical studies has been performed 
to model the energy dependence, or time lags observed in 
the day-scale flares in blazars. One approach considers 
models with a homogeneous emission region (Mastichiadis 
& Kirk 1997; Chiaberge & Ghisellini 1999; Li & Kusunosc 
2000). For instance, Mastichiadis & Kirk (1997) have sug- 
gested that a sudden increase in the maximum energy of 
the accelerated electrons could result in the X-ray and TcV 
flare observed in Mrk 421. Similar results were shown by 
Kataoka et al. (2000), in reproducing the soft-lag observed 
in PKS 2155-304. However, most of these models assume 
a change in some parameter in the emission region (such 
as the magnetic field, maximum energy of electrons, num- 
ber density of electrons, etc.). The observational results 
for Mrk 421 described above lead us to suggest a different 
scenario. 

An alternative scenario has each flare forming as a re- 
sult of a separate electron distribution. This was suggested 
in modeling the high state SED of Mrk 501 (Kataoka et 
al. 1999). They applied the one-zone homogeneous syn- 
chrotron sclf-Compton model to the SED from April 1997, 
and concluded that a single electron distribution is insuf- 
ficient to reproduce the observed synchrotron spectrum. 
This scenario could work for instance when the jet is emit- 
ted intermittently from the central engine, and several sep- 
arate regions generate different emission components. An- 



other viable scenario can be provided by an internal shock 
model, where the light curve results from a superposition 
of many flares due to the collisions of shells, which may 
occur when a faster shell catches up to a slower shell (Ghis- 
ellini 2001, Spada et al. 2001; Sikora et al. 2001) 

In fact, Tanihata et al. (2003) has shown from simula- 
tions that the internal shock model can naturally explain 
various variability properties observed in TeV blazars. 
One observational fact is that the X-ray flares always ap- 
pear to lie on top of an underlying offset-like component. 
In the internal shock model, flares can be considered as 
arising from collisions of shells that had the largest differ- 
ence in the initial velocity. All other collisions generate 
smaller amplitude, longer flares, which pile up to generate 
the offset component. In this case, it was shown that the 
flares will tend to have a higher synchrotron peak energy 
as compared with that of the offset component. 

Assuming that flares are due to an emergence of a new 
component with a higher synchrotron peak energy than 
the pre-existing component, the synchrotron peak energy 
in the observed spectrum would appear as shifting to a 
higher energy. It is interesting to note that our observa- 
tions arc consistent with this, with the overall trend of 
the synchrotron peak energy being higher during flares. 
Furthermore in this case, the spectral evolution would de- 
pend on the relation of the new and pre-existing spectra. 
This would naturally explain the different spectral evolu- 
tion among different flares. 

Finally, we remark again that due to the data gap be- 
tween the EUVE and ASCA data in our observation, the 
results concerning the spectrum peak are somewhat tenta- 
tive. Because of the rapid variability of the source in short 
time scales, the effective area of the instrument is most im- 
portant. Long- look observations by the new observatories 
such as XMM-Newton should provide data allowing a 
more detailed analysis of the flaring mechanism, leading 
to the dynamics of the accelerated electrons in blazar jets. 
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ciety for Promotion of Science for Young Scientists, and 
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Fig. 1. — The joint fit with a quadratic function to the LECS, MECS, GIS, and PCA spectra during the period when the ASCA, BeppoSAX, 
and RXTE observations overlapped (MJD 50927.03 — 50927.34). The top panel shows the photon spectrum for each instrument. The solid 
line shows the best-fit model. The bottom four panels shows the residuals of the observed data to the model. See Table 1 for the values of 
the fit parameters. 

Fig. 2. — The time evolution of the spectrum fit parameters to the combined EUVE and ASCA spectrum, assuming a quadratic function 
(see text). The top three panels show the peak flux (top), peak energy (second), and the curvature parameter (third; see definition in Table 1). 
The flux is in units of 10~ 12 erg cm~ 2 s -1 , and peak energy is in keV. The bottom panel shows the reduced-% 2 for each of the fits. Three 
out of the 56 spectra showed a poor fit (P(x 2 ) <5%) due to several line or edge-like structures, where the fit is significantly improved after 
excluding the energy ranges that included these structures, shown as the crosses. 

Fig. 3. — The correlation of the synchrotron peak flux and peak energy for Mrk 421, assuming a quadratic function. Results derived from 
the combined EUVE and ASCA spectra during the 1998 campaign arc shown as the filled circles. Results from the BeppoSAX spectra from 
the 1998 campaign and 1997 observation are shown as the open circles and squares, respectively. The fluxes for the BeppoSAX data are 
multiplied by a factor of 1.07 to normalize them to the ASCA fluxes. Each data point corresponds to the value derived from a spectrum 
integrated over 10 ksec. 

Fig. 4. — The time evolution of the correlation of the synchrotron peak flux vF v pf ,^ and peak energy -E pea k f° r Mrk 421 during the 
1998 campaign, plotted separately for each of the flares defined in Figure 2. The initial values in each period are shown as a star, and the 
connecting line shows the time order. This Figure illustrates that the relation of J'-F^peak and -Epcak differ from flare to flare. 

Fig. 5. — The spectral evolution during the rise of the second flare in the Mrk 421 observation during the 1998 campaign. The ratio of 
the spectrum to that of the onset of the rise is shown in the lower panel. It can be seen that the flare starts with a hardening at the highest 
energy range with a concave shape in the ratio, which with time, propagates down to the lower energies. In the last spectrum, which is 40 ks 
after the initial one, the whole spectrum is described with the new slope. 

Fig. 6. — The spectral evolution during the rise of the fourth flare in the Mrk 421 observation during the 1998 campaign. The ratio of the 
spectrum to that of the onset of the rise is shown in the lower panel. 

Fig. 7. — The spectral evolution during the rise of the fifth flare in the Mrk 421 observation during the 1998 campaign. The ratio of the 
spectrum to that of the onset of the rise is shown in the lower panel. 

Fig. 8. — The X-ray spectra of Mrk 421 obtained simultaneously with the (a) HEGRA and (b) Whipple observations during the 1998 
campaign. The X— ray spectra with different symbols correspond to each TeV observation shown in the label. The X-ray data include those 
collected by the EUVE and ASCA satellites for the HEGRA pointing, and additionally, by the RXTE satellite for the Whipple observations, 
and clearly show the curvature of the synchrotron spectrum. 

Fig. 9. — The correlation of the synchrotron peak flux and the TeV flux (in Crab units) for the simultaneous observations of Mrk 421 
during the campaign in 1998. The synchrotron peak flux is derived from observations by EUVE and ASCA (and RXTE for the Whipple 
pointings) during the exactly simultaneous epochs where the TeV data are available. The letters H and W in the legend denote HEGRA and 
Whipple pointings. 



Table 1 



Results of Spectral Fits for Mrk 421 during MJD 50927.03 - 50927.34 



Instrument 


Model 




r 2 b 


-Ebr° -Epea,k° 


^,peak d Q g \l (dof) 


LECS (0.12-4) 


PL 


2.20±0.01 






2.85 (83) 


MECS (1.65-10.5) 


PL 


2.76±0.02 






1.28 (193) 


GIS (0.7-7) 


PL 


2.52±0.01 






1.29 (124) 


PCA (3-25) 


PL 


2.91±0.01 






2.63 (50) 


LECS, MECS, GIS, PCA 
LECS, MECS, GIS, PCA 
LECS, MECS, GIS, PCA 
LECS,MECS,GIS,PCA 


broken PL 
curved PL° 
cutoff PL f 
quadratic 5 


2.21±0.01 
1.77±0.04 
2.17±0.01 


2.90±0.01 
1.66±0.03 


2.24±0.04 

2.5±0.3 
7.80±0.14 

0.47±0.02 


1.86 (451) 
1.08 (451) 
1.74 (452) 

217±3 0.42±0.01 1.07 (452) 



a The numbers in parentheses show the energy range used in the spectral fitting 
b Photon index 
c in keV 

d in 10~ 12 crg cm" 2 s" 1 
°F{E) = KE- r i(l + -^-r r2 
1 F(E) = KE~ Tl cxp(-E/E bl ) 

Ho SV F„(E) = log^F^peak) - a(\o S E - log,E pcak ) 2 
Note. — All errors represent the la error. 
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Table 2 

Spectral Fit Parameters for the combined EUVE /ASCA spectra. 



O bscrvation 


E l> 

■^pcak 


J/r f,pcak 


a d 




xl (dof) 


i 


0.53±0.02 


230±2 


0.39±0 


02 


1.06 


(66) 


2 


0.51±0.02 


244±3 


0.39±0 


02 


0.95 


(85) 


3 


0.57±0.02 


266±2 


0.39±0 


01 


0.95 


(122) 


4 


0.69±0.02 


280±2 


0.48±0 


01 


1.22 


(106) 


5 


0.62±0.01 


331±2 


0.42±0 


01 


1.08 


(116) 


6 


0.72±0.02 


372±2 


0.40±0 


01 


0.97 


(137) 


7 


0.63±0.04 


377±5 


0.41±0 


02 


0.89 


(145) 


8 


0.60±0.03 


323±4 


0.46±0 


02 


1.20 


(87) 


9 


0.49±0.02 


277±3 


0.38±0 


02 


0.80 


(82) 


10 


0.55±0.02 


270±2 


0.36±0 


01 


0.92 


(79) 


11 


0.58±0.02 


278±2 


0.36±0 


01 


1.10 


(137) 


12 


0.59±0.02 


307±2 


0.36±0 


01 


1.04 


(126) 


13 


0.59±0.02 


337±2 


0.38±0 


01 


0.97 


(125) 


14 


0.66±0.02 


358±2 


0.39±0 


01 


1.01 


(154) 


15 


0.69±0.02 


324±3 


0.43±0 


02 


0.73 


(97) 


16 


0.65±0.02 


303±2 


0.37±0 


01 


0.89 


(87) 


18 


0.71±0.02 


333±2 


0.35±0 


01 


0.93 


(98) 


19 


0.76±0.02 


364±2 


0.37±0 


01 


0.92 


(173) 


20 


0.68±0.02 


346±2 


0.38±0 


01 


1.16 


(144) 


21 


0.77±0.02 


337±2 


0.42±0 


01 


1.04 


(132) 


22 


0.86±0.02 


391±2 


0.40±0 


01 


1.14 


(169) 


23 


0.89±0.02 


458±2 


0.43±0 


01 


0.97 


(184) 


24 


0.73±0.03 


423±3 


0.34±0 


02 


1.21 


(119) 


25 


0.87±0.03 


379±3 


0.41±0 


02 


0.85 


(113) 


26 


0.85±0.02 


366±2 


0.36±0 


01 


0.88 


(95) 


27 


0.71±0.02 


336±2 


0.36±0 


01 


1.14 


(153) 


28 


0.68±0.02 


397±3 


0.38±0 


01 


1.09 


(146) 


29 


0.74±0.02 


404±3 


0.39±0 


01 


0.85 


(125) 


30 


0.70±0.01 


423±2 


0.41±0 


01 


1.08 


(176) 


31 


0.69±0.01 


442 ±2 


0.40±0 


01 


0.88 


(181) 


32 


0.74±0.02 


439±3 


0.41±0 


01 


0.81 


(124) 


33 


0.66±0.02 


459±3 


0.38±0 


01 


1.07 


(136) 


34 


0.83±0.02 


390±3 


0.43±0 


01 


0.97 


(93) 


35 


0.72±0.02 


350±2 


0.34±0 


01 


1.02 


(135) 


36 


0.73±0.03 


365±2 


0.30±0 


01 


0.78 


(136) 


37 


0.85±0.02 


388±2 


0.34±0 


01 


0.96 


(165) 


38 


0.95±0.02 


413±2 


0.33±0 


01 


0.97 


(187) 


39 


1.11±0.02 


515±2 


0.33±0 


01 


1.07 


(182) 


40 


1.09±0.03 


469±2 


0.36±0 


01 


0.96 


(146) 


41 


1.02±0.04 


423±2 


0.32±0 


02 


1.11 


(153) 


42 


0.88±0.02 


348±2 


0.37±0 


01 


0.70 


(92) 


43 


0.76±0.03 


352±2 


0.32±0 


01 


0.98 


(131) 


44 


0.76±0.03 


369±3 


0.30±0 


02 


1.23 


(84) 


45 


0.94±0.03 


343 ±2 


0.27±0 


01 


0.98 


(97) 


46 


0.80±0.02 


293±1 


0.29±0 


01 


0.96 


(126) 


47 


0.92±0.02 


308±2 


0.33±0 


01 


1.09 


(138) 


48 


0.88±0.03 


311±2 


0.37±0 


01 


1.04 


(130) 


49 


0.85±0.03 


378±3 


0.30±0 


02 


1.13 


(80) 


50 


1.13±0.03 


476±2 


0.32±0 


01 


1.04 


(132) 


51 


1.12±0.02 


443 ±2 


0.35±0 


01 


0.79 


(144) 


52 


0.95±0.03 


416±3 


0.36±0 


02 


0.86 


(112) 


53 


0.95±0.02 


429±2 


0.36±0 


01 


1.04 


(171) 


54 


0.85±0.02 


365±2 


0.38±0 


01 


1.13 


(180) 


55 


0.82±0.02 


341±2 


0.34±0 


01 


1.02 


(152) 


56 


0.98±0.02 


372±2 


0.35±0 


01 


1.10 


(157) 




0.99±0.03 


396±2 


0.32±0 


01 


1.00 


(166) 



a Each number represents the observation divided into 10 ksec seg- 
ments 

b in kcV 

c in 10 _12 crg cm -2 s -1 
d Dcfinition is given in Table 1 



Note. - 



All errors represent the la error. 
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Table 3 

Spectral Fit Parameters for the BeppoSAX spectra 



Observation 41 






Quadratic 






Curved Power-law 








■f'pcak 




a d 


xl (dof) 


^pcak 




Xl (dof) 


1997.4.29 


1 


0.54±0.03 


137±3 


0.48±0.02 


0.99 (126) 


0.56±0.02 


146±6 


0.96 


(125) 




2 


0.46±0.02 


124±2 


0.48±0.02 


1.21 


(126) 


0.49±0.02 


125±4 


1.23 


(125) 




3 


U.oU±U.Ul 


131±2 


0.46±0.01 


1.16 


(126) 


0.53±0.01 


134±3 


1.15 


(125) 


1997.4.30 


1 


0.53±0.03 


132±4 


0.48±0.02 


1.24 


(126) 


0.55±0.04 


127±5 


1.25 


(125) 




2 


U.47±U.0z 


127±3 


0.47±0.02 


0.88 


(126) 


0.50±0.02 


132±5 


0.85 


(125) 




3 


0.55±0.03 


124±3 


0.51±0.02 


0.85 


(126) 


0.58±0.03 


123±4 


0.84 


(125) 




4 


0.52±0.01 


129±2 


0.47±0.01 


0.81 


(126) 


0.55±0.01 


130±2 


0.79 


(125) 


1997.5.1 


1 


0.45±0.03 


129±4 


0.45±0.02 


1.11 


(126) 


0.48±0.03 


131±6 


1 11 


( 125) 




2 


0.57±0.03 


144±3 


0.44±0.02 


1.03 


(126) 


0.59±0.03 


146±4 


1 03 


( 125) 




3 


0.58±0.02 


178±4 


0.49±0.02 


0.84 


(126) 


0.61±0.03 


170±5 


0.85 


(125) 




4 


0.52±0.01 


143±2 


0.44±0.01 


1.03 


(126) 


0.55±0.02 


145±2 


1.01 


(125) 


1997.5.2 


1 


0.57±0.02 


171±3 


0.46±0.01 


0.86 


(126) 


0.60±0.02 


174±4 


0.85 


(125) 


1997.5.3 


1 


0.42±0.03 


120±4 


0.48±0.02 


0.97 


(126) 


0.39±0.06 


116±4 


u.yo 


( 125 ) 




2 


0.42±0.02 


123±3 


0.44±0.02 


1.13 


(126) 


0.41±0.05 


120±4 


1 14 


( 125) 


1997.5.4 


1 


0.31±0.02 


105±4 


0.50±0.03 


1.23 


(126) 


0.22±0.09 


104±5 


1 25 


(125) 




2 


0.32±0.02 


103±3 


0.50±0.02 


1.23 


(126) 


0.32±0.05 


103±4 


1 24 


( 125) 


1997.5.5 


1 


0.38±0.02 


127±4 


0.47±0.02 


1.24 


(126) 


0.42±0.03 


131±5 


1 .24 


(125) 




2 


0.37±0.02 


127±3 


0.46±0.02 


1.24 


(126) 


0.40±0.03 


129±4 


1 .25 


(125) 


1998.4.21- 


1 


0.69±0.04 


341±8 


0.42±0.02 


1.04 


(126) 


0.69±0.04 


357±13 


— i — m — 


(125) 




2 


0.85±0.02 


404±6 


0.42±0.01 


1.12 


(126) 


0.93±0.03 


397±7 


1.06 


(125) 




3 


0.87±0.02 


378±5 


0.37±0.01 


1.03 


(126) 


0.91±0.03 


379±6 


1.04 


(125) 




4 


0.83±0.03 


337±5 


0.37±0.01 


0.93 


(126) 


0.87±0.04 


337±7 


0.94 


(125) 






U. 1 HU.UO 


300±5 


0.34±0.01 


0.90 


(126) 


0.72±0.04 


310±7 


0.86 


(125) 




6 


0.78±0.03 


274±4 


0.37±0.01 


1.31 


(126) 


0.82±0.03 


274±5 


1.33 


(125) 




7 


0.64±0.03 


239±4 


0.33±0.02 


1.04 


(126) 


0.68±0.03 


235±5 


1.02 


(125) 




8 


0.67±0.03 


236±5 


0.38±0.02 


0.95 


(126) 


0.70±0.04 


235±6 


0.96 


(125) 


1998.4.23- 


1 


0.65±0.03 


309±8 


0.44±0.02 


1.11 


(126) 


0.70±0.05 


287±8 


1.06 


(125) 




2 


0.58±0.03 


257±6 


0.38±0.02 


1.07 


(126) 


0.62±0.03 


257±8 


1.08 


(125) 




3 


0.56±0.02 


247±5 


0.40±0.02 


0.90 


(126) 


0.58±0.02 


254±7 


0.88 


(125) 




4 


0.58±0.03 


248±6 


0.43±0.02 


1.04 


(126) 


0.61±0.03 


239±7 


1.01 


(125) 




5 


0.57±0.02 


253±5 


0.44±0.02 


0.86 


(126) 


0.60±0.02 


251±7 


0.87 


(125) 




6 


0.52±0.02 


233±5 


0.44±0.02 


1.17 


(126) 


0.53±0.03 


227±6 


1.18 


(125) 




7 


0.49±0.02 


220±5 


0.43±0.02 


1.34 


(126) 


0.51±0.03 


218±6 


1.36 


(125) 




8 


0.46±0.03 


219±7 


0.38±0.02 


1.07 


(126) 


0.44±0.06 


213±8 


1.08 


(125) 




9 


0.45±0.03 


217±7 


0.42±0.02 


1.04 


(126) 


0.34±0.09 


203±7 


1.01 


(125) 




10 


0.44±0.03 


203±5 


0.40±0.02 


1.14 


(126) 


0.39±0.06 


194±6 


1.14 


(125) 



a Each number represents the observation divided into 10 kscc segments 
b in kcV 

c in 10 — 12 erg cm -2 s _1 

d Definition is given in Table 1 



Note. — All errors represent the la error. 



Table 4 



Simultaneous TeV flux and Synchrotron Peak Energy and Flux for Mrk 421 



Start Time 


Duration 


TeV 






X- 


-ray 




(MJD) 


(kscc) 


Observatory 


TeV flux" 


Sp Ca k b 




a d 


Xl (dof) 


50926.8814 


11.0 


HEGRA 


0.46±0.10 


0.48±0.03 


226±6 


0.45±0.04 


0.75(33) 


50927.1390 


13.7 


Whipple 


0.33±0.09 


0.47±0.01 


218±2 


0.42±0.01 


0.95(216) 


50927.8839 


10.5 


HEGRA 


0.79±0.12 


0.74±0.02 


381±3 


0.43±0.01 


1.18(210) 


50928.2273 


5.3 


Whipple 


0.33±0.07 


0.56±0.02 


269±3 


0.42±0.01 


0.95(142) 


50928.8824 


10.5 


HEGRA 


0.58±0.10 


0.66±0.02 


316±2 


0.41±0.01 


0.90(221) 


50929.8848 


10.1 


HEGRA 


1.10±0.13 


0.77±0.02 


424±2 


0.36±0.01 


1.13(205) 


50930.1516 


12.2 


Whipple 


0.54±0.05 


0.81±0.02 


364±2 


0.40±0.01 


0.99(274) 


50930.8845 


9.8 


HEGRA 


1.41±0.15 


0.66±0.03 


472±4 


0.39±0.01 


1.15(244) 


50931.1467 


10.3 


Whipple 


0.58±0.14 


0.72±0.02 


367±3 


0.38±0.01 


0.94(229) 


50932.2041 


6.9 


Whipple 


0.67±0.07 


0.84±0.02 


367±2 


0.34±0.01 


1.31(239) 


50933.2358 


3.4 


Whipple 


0.67±0.11 


0.90±0.02 


454±2 


0.38±0.01 


1.21(238) 



a in Crab units 
b in kcV 

c in 10- 12 erg cm" 2 s" 1 
d Dcfinition is given in Table 1 



Note. — All errors represent the la error. 



Evolution of spectrum of Mrk 421 
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Evolution of spectrum of Mrk 421 
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Energy [keV] 



• 50926.8 
O 50927.8 
■ 50928.8 

□ 50929.3 

♦ 50930.8 
O 50927.1 
A 50928.2 
A 50930.1- 
▼ 50931.1! 
V 50932.2 

* 50933.2 
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